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ABSTRACT 

We present the results of new AGILE observations of PSR B 1509-58 per- 
formed over a period of ~2.5 years following the detection obtained with a subset 
of the present data. The modulation significance of the lightcurve above 30 MeV 
is at a 5a confidence level and the lightcurve is similar to those found earlier by 
COMPTEL up to 30 MeV: a broad asymmetric first peak reaching its maximum 
0.39 ± 0.02 cycles after the radio peak plus a second peak at 0.94 ± 0.03. The 
gamma-ray spectral energy distribution of the pulsed flux detected by COMP- 
TEL and AGILE is well described by a power-law (photon index a = 1.87±0.09) 
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with a remarkable cutoff at E c = 81 ±20 MeV, representing the softest spectrum 
observed among gamma-ray pulsars so far. The pulsar luminosity at E > 1 MeV 
is L 7 = 4.21q 2 x 10 35 er §/ s ; assuming a distance of 5.2 kpc, which implies a spin- 
down conversion efficiency to gamma-rays of ~ 0.03. The unusual soft break in 
the spectrum of PSR B 1509-58 has been interpreted in the framework of polar 
cap models cLS cL SI gnature of the exotic photon splitting process in the strong 
magnetic field of this pulsar. In this interpretation our spectrum constrains the 
magnetic altitude of the emission point (s) at 3 km above the neutron star surface, 
implying that the attenuation may not be as strong as formerly suggested because 
pair production can substitute photon splitting in regions of the magnetosphere 
where the magnetic field becomes too low to sustain photon splitting. In the 
case of an outer-gap scenario, or the two pole caustic model, better constraints 
on the geometry of the emission would be needed from the radio band in order 
to establish whether the conditions required by the models to reproduce AGILE 
lightcurves and spectra match the polarization measurements. 

Subject headings: gamma rays: observations — pulsars: general — pulsars: in- 
dividual (PSR J1513-5908 (B1509-58)) — stars: neutron 



Introduction 



PSR B1509-58 (J1513-5908) was discovered as an X-ray pulsar with t he Einstein satellite 
during an observation of the supernova remnant (SNR) MSH 1 5-52 (jSeward fc Harnderi 
19821 ) . The source was soon also detected at radio frequencies by Manchester et al.l (119821 ) , 
with a derived distance suppo rting the association with the SNR (d ~ 5.2 kpc, as calculated 
from HI measurements from Gaensler et al. 1999 and in agreement with the most recent 
distance derived using the dispersion measurelfwith a period P H 150 ms and a period 
derivative P ~ 1.53 x 10~ 12 s s^ 1 , assuming the standard dipole vacuum model, the estimated 
spin-down age for this pulsar is 1570 years (among the shortest for radio pulsars) and its 
inferred surface magnetic field is one of the highest observed for an ordinary radio pulsar: 
B = 3.1 x 10 13 G, as calculated at the pole B Its rotational energy loss rate is E = 
1.8 x 10 37 erg/s. 



1 see the ATNF Pulsar Catalogue ( jhttp: / / www. atnf . csiro . au / research / pulsar /psrcat /) for the updated 
distance measurement derived from the dispersion measure. 

2 The magnetic field strength at the pole is twice the value quoted in the ATNF Pulsar Catalogue. 
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PSR B 1509-58 and its nebula have been extensively observed in the X-ray energies 
since the Eighties with the Einstein and EXOSAT satellites . The detect i on of pulsed emis- 
sion in the hard X-rays dates back to the early Nineties ( Kawai et al. 199lh with Ginqa 
in th e 2-60 keV energy range. During a 20 yrs-long radio monitoring ((Livingstone et aL 



20051 ). PSR B1509-58 has not shown any glitch activity, at variance with the general be- 



havior of young radio p ulsars, which usually show some glitch activity. The analysis of 
Livingstone et al.l ( 120051 ). using radio and X-rays (collected with RossiXTE), yielded a very 
accurate measurement of the braking index of n = 2.839 ± 0.003, close to the canonical 
value n = 3 for braking by magn etic dipo le radiation in vacu um alone. Obser yations with 



the R OSAT jTrussoni et al.lll996h . ASCA jTamura et alill996h and BeppoSAX jMineo et al. 



20011 ) satellites were performed in the Nineties, characterizing the spectrum of the pulsed 
emission and the morphology of the remnant as possibly due to the presence of several 
components, interacting via collimated outfl ows from the pulsar. The nebula has been ex- 
tensively observed with the Chandra satellite ( iGaensler et al.ll2002l) and it s emission has been 
found up to the TeV energies, with CANGAROO first (ISako et al.ll2000l ) and more recently 
by H.E.S.S. (Uharonian et al-lbood ). 



The young age and the high rotational energy loss rate made this pulsar a promis- 
ing target for the first generation of gamma-ray satellites. In fact, the instruments on the 
Compton Gamma-Ray Observatory (C GRO) observed its pu lsation at lo w gamma-ray ener- 



gies: up to E ~ 7 00 keV with BATSE (jWilson fc et al.lll993h and OSSE (Ul mer et al 



1993 



Matz et al|ll99J ), and in the 0.75 - 30 MeV band with COMPTEL (IKuiper et al.lll999h . 
but it was not detected with high significance by the Energetic Gamma-Ray Experiment 
Telescope (EGRET), the instrument operating at the energies from 30 MeV to 30 GeV. This 
was remar kable, since all ot h er known gamma-ray pulsars show spectral turnovers well above 
100 MeV ([Thompson! 12004 ) . lHarding et al.l (119971 ) suggested that the break in the spectrum 
could be inter preted as due to in hibition of the pair-production caused by the photon-splitting 
phenomenon (lAdler et al.lll970l ). The photon splitting appears, in the frame of the polar cap 
mo dels, in relation with a very high magnetic field. An alternative explanation is proposed 
by IZhang fc Cheng! (120001 ) using a three dimensional outer gap model. They propose that 
the gamma-ray emission is produced by synchrotron-self Compton radiation above the outer 
gap. 

Ten years after the CGRO, the observation in the gamma-ray band are possible again 
with the advent of two high-energy missions: the Astro-rivelatore Gamma a Immagini 
LEggero (AGILE) and Fermi satellites. With its large field of view (~ 2.5 sr) and its 
< 1/xs time tagging cap ability (jTavani et al.ll2009i ) (with 200/xs absolute timing accuracy, see 



Pellizzoni et al.l l2009bl ) . AGILE is very well suited for the observation of pulsars between 



30 MeV and a few GeV with its Gamma-Ray Imaging Detector (GRID). In particular, de- 
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spite its lower sensitivity in the GeV band, the GRID on board AGILE has an effective area 
below 100 MeV (~200 cm 2 at 50 MeV) comparable with that of Fermi. AGILE obtained 
the first detection of PSR B1509-58 in the EGRET band JPellizzoni et al.ll2009ah confirm- 
ing the occurrence of a spectral break, although no precise flux measurements were possible 
due to low counts statistics. Recently Fermi also reported its detection of PSR B1509-58 
( lAbdo et al.ll2010ai ). In this paper we present the results of a ~ 2.5 yr monitoring campaign 
of PSR B 1509-58 with AGILE, that improved counts statistics, and therefore lightcurve 
characterization, with respect to earlier AGILE observations. With these observations the 
spectral energy distribution (SED) at energies E < 300 MeV, where the remarkable spectral 
turnover is observed, can be assessed. 



2. AGILE Observations, Data Analysis and Results 

PSR Bl 509-58 is within the same region of the sky as the Vela pulsar, an area to which 
AGILE devoted a large amount of observing time (for d etails on AGILE obsery ing strategy, 
timing calibration and gamma-ray pulsars analysis see IPellizzoni et~aUl2009bl laT). Gamma- 



ray photons for this pulsar were collected and analyzed starting in July 2007, up to late 
October 2009 when AGILE started observing in spinning mode due to reaction wheel failure 
0. The large AGILE effective area and long observing time (~ 260 days on target) provide 
a total exposure of 3.8 x 10 9 cm 2 s (E > 100 MeV) during this 2.5 yr period which gives our 
observations a good photon harvest from this pulsar. 



2.1. Timing Analysis 

Simultaneous radio observations of PSR B1509-58 with the 64 m Parkes radiotelescope 
in Australia are ongoing since the epoch of AGILE' s l aunch (MJD 54220; 2007 April 30), 



as part of a timing project for the gamma-ray satellites (IWeltevrede et al.ll2010bl ). and cover 
all of AGILE 's observations. A total of 47 pulsar time of arrivals (ToAs) were collected 
between April 2007 (MJD 54220) and February 2010 (MJD 55233), leading to a r.m.s. of 
the residuals of 900 /is, showing the goodness of the timing model that allowed accurate 
pulse phase tagging of the gamma-ray photons. No glitch was detected in the radio analysis. 
Strong timing noise was present, as expected from a young pulsar, and it was accounted 
for using the fitwaves technique developed in the framework of the TEMP02 radio pul- 



3 This failure is not affecting AGILE/GRID sensitivity and pulsar observations although the new spinning 
mode required calibration revisions that are still ongoing. 
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sar timing software (IHobbs et al.ll2004j . 120061 ) . Using the radio ephemeris provided by the 
Parkes tel escope, we performed th e folding of the gamma-ray lightcurve including the wave 



terms (see iPellizzoni et al.ll2009bl ). An optimized analysis followed, aimed at cross-checking 
and maximization of the significance of the detection, including an energy-dependent event 
extraction angle around the source position based on the instrument point-spread-function 
(PSF). Only high confidence gamma-ray photons (G) were used for the timing analysis of 
this pulsar. The chi-squared (x 2 )-test applied to the 10 bin lightcurve at E > 30 MeV gave 
a detection significance of a = 4.8. The unbinned Z^-test applied to the photons' arrival 
times gave a significance of a = 5.0 with n = 2 harmonics. The difference between the radio 
and gamma-ray ephemerides was AP ra ^ 0j7 < 10 -9 s, well within the error in the parameters, 
showing perfect agreement among radio and gamma-ray ephemerides as expected, further 
supporting our detection and AGILE timing calibration. 

We observed PSR B 1509-58 in three energy bands. We obtained 1210 ± 400 pulsed 
counts (~ 5% of the total source, diffuse gamma-ray emission and residual particle back- 
ground counts) at energies 30 < E < 100 MeV, 820 ± 360 pulsed counts (~ 7% of the total 
counts) at energies 100 < E < 500 MeV. The pu lsed flux was computed considering all the 
counts above the minimum of the lightcurve (see IPellizzoni et al.ll2009bl ) We did not detect 
pulsed emission at a significance a > 2 for E > 500 MeV and thus we ca n only give an upper 
limit at la of < 270 pulsed counts. This is consistent with the fact that lAbdo et al.l ( I2010af ) 
report only a 1.4<r detection at 0.3 < E < 1 GeV with the Fermi data. 

The gamma-ray lightcurves of PSR B1509-58 for different energy bands are shown in 
Fig. [TJ The AGILE lightcurve above 30 MeV shows two peaks at phases 0i = 0.39 ± 0.02 
and 02 = 0.94 ± 0.03 with respect to the single radio peak, here put at phase 0, as obtained 
from the Parkes ephemeris. The peak positions and widths in term of phase are calculated 
using a Gaussian fit, yielding a FWHM of 0.29(6) for the first peak and of 0.13(7) for the 
second peak, where we quote in parentheses (here and throughout the paper) the la error 
on the last digit. The errors considered are statistical, as the systematic errors do not affect 
the measurements of the pulsed counts. The first peak is coincident in phase with the X-ray 
single broad peak and with COMPTEL peak (see lKuiper et al.lll999l and references therein). 
In its highest energy band (10-30 MeV) COMPTEL showed the hint of a second peak (even 
though the modulation had low significance, 2.1cr), which is also visible in the lightcurve 
derived from EGRET observations (30-100 MeV), despit e the fact th at EGRET did not 
have a significant detection of the pulsar above 100 MeV (IFierrol Il995l ) . This second peak 
is coincident in phase with AGILE'S second peak (Fig. [T]) in its lower energy band while it 
appears slightly shifted at energies above 100 MeV. A possible explanation for this shift is 
discussed in section § 3. AGILE thus confirms the previously marginal detection of a second 
peak, statistically significant at 5a, calculated using a chi squared statistics test. 
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2.2. Spectral Analysis, Spectral Energy Distribution 

Based on our exposure, calculated by the GRID scientific analysis task AG_ExmapGen, 
we derived the gamma-ray flux from the number of pulsed counts. This method, though 
typically giving higher statistical errors than the likelihood analysis, is more accurate and 
sensitive to evaluate the flux of this pulsar, given its soft spectrum (and the correspondingly 
large PSF) and the contribution from other nearby and brighter sources and possibly from 
the pulsar wind nebula (PWN), that all affect the spatial analysis. Using this method, the 
problem of modeling the background is dealt with by discarding the counts below the pulsed 
threshold, so that the observed pulsed counts belong to the pulsar. All the other sources 
of systematical errors, such as the effective area and AGILE 's PSF, as well as the diffuse 
emission, contribute to < 10%, so that they are much lower than the statistical errors and 
they are not quoted in the flux measurements. We divided AGILE bandwidth into three 
energy intervals: 30-100 MeV, 100-500 MeV and above 500 MeV. The pulsed fluxes thus 
obtained were F 1 = 10(3) x 10~ 7 ph cm" 2 s" 1 in the 30-100 MeV band, F 7 = 1.8(8) x 
10~ 7 ph cm -2 s -1 in the 100-500 MeV band and a la upper limit F 7 < 8x 10~ 8 ph cm" 2 s~ x 
for E > 500 MeV. Finally, from the total number of pulsed counts we obtained a pulsed flux 
at E > 30 MeV F 7 = 12(2) x 10" 7 ph cm" 2 s" 1 for E > 30 MeV. 

Fig. 2 shows the SED of PSR Bl 509-58 based on AGILE'S and COMPTERS observed 
fluxes. Fermi upper limits are also shown, which are consistent with our measurements at 
a 2a confidence level. COMPTEL observed this pulsar in three energy bands: 0.75-3 MeV, 
3-10 MeV, 10-30 MeV, suggesting a spectral break between 10 and 30 MeV. AGILE pulsed 
flux confirms the presence of a soft spectral break. As shown in Fig. [21 we modeled the 
observed COMPTEL a nd AGILE fluxes wi th a power-law plus cutoff fit using the Minuit 



minimization package (I James &: Rooslll975l ): F(E) = k x E~ a exp[— (E/E C )P], with three 
free parameters: the normalization k, the spectral index a, the cutoff energy E c and allowing 
(3 to assume values of 1 and 2 (indicating either an exponential or a superexponential cutoff). 
No acceptable x 2 values were obtained for a superexponential cutoff, the presence of which 
can be excluded at a 3.5a confidence level, while for an exponential cutoff we found x 2 = 3-2 
for v = 2 degrees of freedom, corresponding to a null hypothesis probability of 0.05. The 
best values thus obtained for the parameters of the fit were: k = 1.0(2) x 10 -4 s _1 cm -2 , 
a = 1.87(9), E c = 81(20) MeV. 

We performed an analysis of the ratio between the two peak heights. The second 
peak appears in the COMPTEL band 10-30 MeV and is observed with AGILE up to E < 
500 MeV: it is harder than the first peak in the COMPTEL energy band, and it is present at 
all energies in the AGILE energy band, so that it might possibly be harder even at AGILE'S 
energies but the low statistics at high energies do not allow us to discriminate. 
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As a consistency check for the pulsed fluxes reported above, a maximum likelihood anal- 
ysis in a region of 10 degrees around the source position was performed to assess possible 
unpulsed contribution from the PWN0, although no detection was reported in t he First Cat- 



alog of High- Confidence Gamma-r ay Sources detected by the AGILE satellite ( IPittori et al. 



20091 ). The likelihood analysis (see iMattox et al.lll996l and for Agile in particular, details will 



be provided in Chen et al. in preparation) took into account the numerous sources present in 
this crowded region (inc luding the extremely bright nearby gamma-ray pulsar J1509-5850, 



Weltevrede et al.ll2010al ). The upper limit found in the AGILE energy range by likelihood 
analysis (F 1 < 40 x 10~ 8 ph cm -2 s _1 above 100 MeV) is above the corresponding pulsed 
flux above 100 MeV (F 7 = 21(6) x 10~ 8 ph cm -2 s^ 1 ). This is compatible with the fact 
that the timing analisys is expected to have for this target a better sensitivity (with respect 
to the likelihood analysis). It is worth noting that PSR B 1509-58 is also not detected by 
likelihood analysis by Fermi (lAbdo et al.ll2010al ) apart from t he > 1 GeV energy ban d where 
the emission could be related to the PWN seen by H.E.S.S. (lAharonian et al.ll2005l ). 



Discussion 



Pulsar magnetosphere models are usually divided into t wo categories, depending on 
the sites for the high-energy emission. In polar cap models (iDaugherty &; Harding! Il996l ) 
the emission comes from the reg ions near the neutron star surface, while outer gap models 
(ICheng et al.lll986l ; lRomanilll996l ) predict that the emission be originated in the regions close 
to the light cylinder. Alternative models predict an emission zone encompassing the whole 
magnetosph ere, which departs from the external rim of the polar cap region: these are the slot 
gap models (jMuslimov fc Harding! 120031 ); o thers predict em ission from alternative regions in 
the outer magnetosphere, the annular gap (IDu et al.ll2010l and references therein). Different 
models predict different spectral and geometrical properties. The bulk of the spin-powered 
puls ar flux is usually em itted in the MeV-GeV en ergy band with spectral breaks at < 10 GeV 
fsee lAbdo elTall boiObl or e.g. lAliu et alJliooi ). PSR B1509-58 has the softest spectrum 
observed among gamma-ray pulsars, with a sub-GeV cutoff at E — 0.08(2) GeV. The second 
softest spectrum and lowest energy cutoff (0 .7(5) GeV) is that of PSR B0656+14, recently 
observed by Fermi (IWeltevrede et al.l l2010al ). The observed lightcurve of PSR B 1509-58 
shows two peaks lagging the radio peak by, respectively, 4>i — 0.39±0.02 and 02 = 0.94±0.03. 
In the following we discuss how the new AGILE observations can constrain the models for 
emission from the pulsar magnetosphere. 



4 See e.g. Pellizzoni et al. ( 2O10l) as an example of study of gamma-ray PWN with AGILE. 
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When PSR B1509-58 was detected in soft gamma-rays but not significantly at E > 
30 MeV, it was proposed t hat the mechanism r esponsible for this low-e nergy spectral bre ak 



might be photon splitting ( [Harding et al.lll997l ). The photon splitting (lAdler et al.lll970f ) is 



an exotic third-order quantum electro-dynamics (QED) process expected when the magnetic 
field approaches or exceeds the critical value defined as B cr = mlc 3 /(eh) = 4.413 x 10 13 G, 
above which quantum effects become relevant. Most current theories for the generation of 
coherent radio emission in pulsar magnetospheres require formation of an electron-positron 
pair plasma developing via electromagnetic cascades. In very high magnetic fields the forma- 
tion of pair cascades can be altered by the process of photon splitting: 7 — > 77, which will 
operate as an attenuation mechanism in the high-field regions near pulsar polar caps. Since 
it has no energy threshold, photon splitting can attenuate photons below the threshold for 
pair production, thus determining a spectral cutoff at lower energies. This process cannot 
operate in the low fields of outer gap models because it only has appreciable reaction rates 
when the magnetic field is at least a significant fraction of the quantum critical field B cr 
(the attenuation coefficient T sp scaling as T sp cx (B/B cr ) 6 = B m ), and magnetic fields strong 
enough are not present in the outer magnetosphere as B ~ r~ 3 . 

In the case of PSR Bl 509-58 a polar cap model with photon splitting would be able 
to explain the soft gamma-ray emission and the low energy spectral cutoff, now quantified 
by AGILE observations. Since the mechanism of photon splitting is, as stated, strongly 
dependent on the magnetic field strength, if the field strength at the emitting region is 
B' > 0.3 (i.e. at heights below 1.3 neutron star radii, Rns), the photon splitting is the 



domi nant means of attenuation that inhibits efficient pair cascade production (IHarding et al. 



19971 ) and then gamma-ray emission. Based on the observed cutoffs, which are related to the 
photons' saturation escape energy, we can derive constraints on the magnetic field strength 
at emission, in the framework of photon splitting: 



f sat ~ 

c esc 



0.077(5' sin 6 kB ,i 



-6/5 



where e s e f c is the photon saturation escape energy and 9kB,o is the angle between the 
photon momentum and the magn etic field vectors at the surface and is here assumed to 
be very small: 9kB,o ^ 0.57° (see IHarding et al.lll997| ). Using the observed energy cutoff 
(e s e f c ~ E = 80 MeV) we find that B' > 0.3, which implies an emission altitude < 1.3Rns, 
which is the height where possibly also pair prod uction could ensue. Th i s alti tude of emis- 
sion agrees with th e polar cap m o dels (see e.g. iDaugherty &: Harding! Il996[ ). A smaller 
energy cutoff, as in IHarding et al.l ( 119971 ). would have implied even lower emission altitude 
and a sharper break, possibly caused by the total absence of pair production. It is appar- 
ent that small differences in the emission position will cause strong differences in spectral 
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shape. This is possibly the reason for the different emission properties of the two peaks as 
observed in the total (AGILE plus COMPTEL) gamma-ray energy band. Also, a trend can 
be observed, from lo wer to higher energi es (see the X-ray lightcurve for the trend in the first 
peak, as in Fig. 3 of iKuiper et al.lll999l ). of the peaks slightly drifting away from the radio 
peak. This we assume to be another signature of the fact that small variations in emission 
height can be responsible for s ensible changes in th e lightcurves in such a high magnetic 
field. The scenario proposed by lHarding et al.l (119971 ) is strengthened by its prediction that 
PSR B0656+14 should have a cutoff with an intermediate value between PSR B1509-58 and 
the other gamma-ray pulsars. The main reason for the parallel between the two pulsars was 
at the time the fact that they had, respectively, the highest and second highest inferred mag- 
netic fields. At present, however, there are a handful of gamma-ray pulsa rs with higher ma; 



netic field than PSR B0656+14 in the Fermi First Year Pulsar Catalog (A bdo et al.ll20 10b 



which do not show a low energy cutoff. None theless, PSR B1509- 5 8 (see 



Crawford et all 120011 ) and PSR B0656+14 foe Luca et all 120051 : IWeltevrede et al 



Kuiper et a 



.111999 : 
2010ah 



both show evi dence of being aligned ro tators, which could imply polar cap emission, as is 
also hinted by iBai fc Spitkovskyi (120101 ) . 



A soft cutoff (below 1 GeV) is in principle possible for polar cap scenarios even without 
invoking photon splitting attenuation. In polar cap models the strong magnetic field permits 
one-photon pair creation that attenuates super-GeV photons in Crab-like (e.g. PSR B1509- 
58, based on the parameter B / P 2 ) and Vela-like pulsars (e.g. PSR B0656+14), whereas pair 
creation in outer gap models is mediated through the two-ph oton process involving surface 



thermal X-rays as targets. According to the calculations of iDaugherty fc Harding! ( 119961 ). 
emission from the regions close to the polar caps is possible when a ~ 6b, where a is the 
angle between the rotation and the magnetic axis and 8b is the half-angle of the gamma- 
beam emerging from the polar cap. Furthermore, with emission from the polar caps, or some 
(> 2) polar cap radii, the pulse profile at high energies can have either one (as in the case 
of PSR B0656+14) or two peaks, with a peak-to-peak phase separation as large as 0.4-0.5 
(albeit slightly smaller than what is observed for PSR B1509-58 at the highest energies). 

The polar cap model as an emission mechanism is nowadays debated. On one hand 
theoretical objections ar ise from the fact that t he angular mom entum is not conserved in 
polar cap emission (see ICohen fc Treves! 1 19721 . iHollowayl 119771 . Treves et al. in prepara- 
tion). At the same time, mounting evidence of a preferential explanation of the observed 
gamma-ray lightcurves with hi gh altitude cascade s is also coming from the recent results 
by the Fermi satellite (see e.g. lAbdo et al.ll2010bl ). In the case of PSR B1509-58, the de- 
rived gamma-ray luminosity from the flux at E > 1 MeV, considering a 1 sr beam sweep is 
L 7 = 4. 21^2^5 2 x 10 35 er g/ s ) where d^,2 indicates the distance in units of 5.2 kpc. While 
traditionally the beaming fraction (/n) was considered to be the equivalent of a 1 sr sweep, 
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nowadays (see e.g. IWatters et al.ll2009[ ) the tendency is to consider a larger beaming fraction 
(fn ~ I); close to a 4ti sr beam. Using /q = 1 in our calculations, we would have obtained 
L 7 = 5.8ig'g(i5 2 x 10 36 erg s -1 . Thus the maximum conversion efficiency of the rotational 
energy loss (E « 1.8 x 10 37 erg s _1 , see § 1) into gamma-ray luminosity is 0.3. Our result is 
not easily comparable with the typical gamma-ray luminosities above 100 MeV, because for 
PSR Bl 509-58 this energy band is beyond the spectral break. Using AGILE data alone we 
obtained a luminosity above 30 MeV L 7 = 5.2(6)d§ 2 X 10 35 erg/s, again for a 1 sr beam. If 
the gamma-ray luminosity cannot account for a large frac tion of the rotational en ergy loss, 
then the angular momentum conservation objection from I Cohen fc Treves! (119721 ) becomes 
less cogent for this pulsar, exactly as it happens for the radio emission. For PSR B0656+14 
no outer magnetosphere emission model seemed to satisfy the observed features and a lower 
magnetosphere model, with an aligned geometry between the rotational and magnetic axes, 
has been proposed and seems plausible from polarization studies. Its efficiency in the con- 
version of the rotational energy l oss into gamma-ray luminosity is one of the lowest observed 
for the gamma-ray pulsars (see iPellizzoni et al.l l2009al lAbdo et al.l l2010bl ): rj = 0.01, not 
violating the constraints imposed by the conservation of angular momentum. 

Alternatively, if such an efficiency as that of PSR B1509-58 were incompatible with this 
conservation law, an interpretation of PSR B1509-58 emission should be so ught in the frame 
of the three dimensional outer magnetosphere gap model, as was done by IZhang fc Cheng 
(120001 ). According to their model, hard X-rays and low energy gamma-rays have the same 
origin: they are produced by synchrotron self-Compton radiation of secondary electron- 
positron pairs of the outer gap. Therefore, as observed, the phase offset of hard X-rays and 
low energy gamma-rays with respect to the radio pulse is the same, with the possibil it y of a 
small lag due to the thickness of the emission region. According to IZhang fc Cheng! (120001 ) 
estimates, a magnetic inclination angle a ~ 60° and a viewing an gle ( ~ 75° ar e required to 



reproduce the observed lightcurve. Similarly, for PSR B0656+14, IWeltevrede et al.l ( 12010a! ) 
argue that large a and ( angles are required to reproduce the observed light c urve in the 



framework of outer gap models. Finally, using the simulations of IWatters et al.l (120091 ). who 
produced a map of pulse profiles for different combinations of angles a and ( in the different 
emission models, the observed lightcurve from AGILE i s best reproduced if a ~ 35° and 
( ~ 90°, in the framework of the two pole caustic model (IDyks fc Rudakil2003l ). 



Since the parameters used for the application of the outer gap model to PSR B 1509-58 
were based on its former observations by COMPTEL, AGILE sp ectrum does not precisely 
fit the spectrum predicted by the model of IZhang fc Cheng! (120001 ) . Furthermore, the values 
of a and ( required by this model are not in good agreemen t with the corresponding values 
obtained with radio measurements. In fact, ICrawford et al.l (120011 ) observe t hat a must be 



< 60° at the 3a confidence level. The prediction obtained by the simulations of IWatters et al. 
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( 120091 ) for slot gap emission is in better agreement with the radio polarization observations 
than what predicted in the outer gap framework. In fact, in the framework of the rotating 



vector model (RVM, see e.g. lLorimer fc Kramerll2004l and references ther ein uCrawford et al. 



( I200ll ) also propose that, if the restriction is imposed that £ > 70° (IMelatosI 119971 ) then 
a > 30° at the 3a level. For these values, however, the Melatos model for the spin down of 
an oblique rotator predicts a braking index n > 2.86, slightly inco nsistent with the observed 
value (n = 2.839(3), see §1). Also in the case of PSR B0656+14, IWeltevrede et all (boiOah 
conclude that the large values of a and ( are somewhat at odds with the constraints from 
the modeling of the radio data and the thermal X-rays which seem to imply a more aligned 
geometry. Improved radio polarization measurements would help placing better constraints 
on the pulsar geometry and therefore on the possibility of a gap in the extended or outer 
magnetosphere, but the quality of the polarization measurements from lCrawford et al.l (1200 ll ) 
is already excellent, the problem being that PSR B1509-58, like most pulsars, only shows 
emission over a limited pule phase range and therefore the RVM models are highly degenerate. 
At present the geometry privileged by the state of the art measurements is best compatible 
with polar cap models. Higher statistics in the number of observed gamma-ray pulsars could 
help characterize a class of "outliers" having gamma-ray emission from the polar caps, which 
potentially constitute a privileged target for AGILE. 



4. Conclusions 



In this paper we present the result of a 2.5 yr long observation campaign of PSR B1509- 
58 with AGILE. With respect to our previous work (jPellizzoni et al.ll2009ai ) the increased 
statistics allowed us to perform an improved lightcurve analysis and to better constrain the 
soft spectral cutoff observed for this pulsar. 



i) Using the Parkes radio ephemeris, AGILE firmly confirmed the detection of gamma-ray 
pulsation with good significance (~ 5a) from PSR B1509-58. 

ii) The observed lightcurve shows two peaks which lag the radio peak of, respectively, 
0.39 ± 0.02 and 0.94 ± 0.03 cycles, as obtained from a Gaussian fit of the peaks. 
PSR B1509-58 presents a single peaked profile up to energies E > 10 MeV where 
COMPTEL detected an additional peak with lower significance. AGILE confirmed 
the existence of this second harder peak in the 30-500 MeV energy band. 

iii) The detection of pulsed emission by AGILE at E > 30 MeV, confirming the presence 
of a soft spectral break, moves the cutoff slightly up, to E ~ 80 MeV, in agreement 
with the previous estimates of a cutoff at energies below 100 MeV. 
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Our observations are compatible with emission from the polar cap regions powered by photon 
splitting cascades. This interpretation could represent a physical measurement related to the 
QED photon splitting process. The fact that polar cap emission at high energies appears 
rare might be explained by the requirement that a number of conditions concur to have low 
magnetosphere emission: an aligned geometry and a high magnetic field, without conflicting 
with the conservation laws. With the AGILE capability of observing with good sensitivity 
at E > 30 MeV, it will be interesting to investigate the highly magnetized pulsars population 
as a possible contributor to a new class of "soft" gamma-ray pulsars. Alternative emission 
models rely on better knowledge of the geometry of PSR B1509-58. 
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Fig. 1. — Phase-aligned gamma-ray light-curves of PSR B1509-58. Radio main peak is at 
phase 0. The start of the y-axis coincides with the minimum of the pulsed fraction and, 
consequently, with the background level. From top to bottom: AGILE high energy band 
(> 100 MeV), 20 bins, 7.5 ms resolution; AGILE "soft" energy band (< 100 MeV), 10 
bins, 15 ms resolution; CO MP TEL high energy ba nd (10-30 MeV) and COMPTEL whole 
bandwith (0.75-30 MeV) ffrom lKuiper et alill999h . 
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Fig. 2. — Spectral energy distribution of PSR B1509-58 (solid line) obtained from a fit of 
pulsed fluxes from soft to hard g amma rays. The three round points represent COMPTEL 
observations (IKuiper et elI.I Il999f ) . The two square points represent AGILE pulsed flux in 
two bands (30 < E < 100 MeV and 100 < E < 500 MeV). The red horizontal bar and 
arrow emerging from it repres ent A GILE upper l imit above 500 MeV. The two green arrows 
represent Fermi upper limits (lAbdo et al J 1201 Oaf ) 



